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Ferroelectricity in multiferroic magnetite Fe3C>4 driven by noncentrosymmetric 
Fe 2+ /Fe 3+ charge-ordering: First-principles study 

Kunihiko Yamauchi 1 , Tetsuya Fukushima 1 , and Silvia PicozziQ 

1. Consiglio Nazionale delle Ricerche - Istituto Nazionale di Fisica della Materia (CNR-INFM), 
CASTI Regional Lab., 67100 L'Aquila, Italy 

(Dated: June 2, 2009) 

By means of first-principles simulations, we unambiguously show that improper ferroelectricity 
in magnetite in the low-temperature insulating phase is driven by charge-ordering. An accurate 
comparison between monoclinic ferroelectric Cc and paraelectric P2/c structures shows that the 
polarization arises because of "shifts" of electronic charge between octahedral Fe sites, leading to a 
non-centrosymmetric Fe 2+ /Fe 3+ charge-ordered pattern. Our predicted values for polarization, in 
good agreement with available experimental values, are discussed in terms of point-charge dipoles 
located on selected Fe tetrahedra, pointing to a manifest example of electronic ferroelectricity driven 
by charge rearrangement. 

PACS numbers: Valid PACS appear here 
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"Improper multiferroics" p], 0, 0, 0, Q are attractive 
multifunctional materials, where magnetism and ferro- 
electricity are strongly coupled. They can be classi- 
fied according to different driving forces, which primar- 
ily break the spatial inversion symmetry paving the 
way to ferroelectricity : 1) spin order, 2) charge order 
(CO) and 3) orbital order. In this letter, we will fo- 
cus on the second group, where LuFe204 has emerged 
as a prototype. @ However, recently Fes04 has been 
also suggested as a ferrimagnet with ferroelectricity be- 
ing induced by charge-ordering, [7j which would depict 
magnetite as the first multiferroic known to mankind. 
FesC^shows the well-known first order metal-insulator 
(Verwey) transition at Ty ~120K, below which the crys- 
tal structure changes from cubic Fd3m to monoclinic 
structure and Fe 2+ /Fe 3+ charge ordering is observed at 
the Fe B sites in the inverse-spinel AB2O4 lattice. 0, Q 
Earlier experiments [Io[ suggested Fes04 to show a spon- 
taneous polarization at low temperatures, with the struc- 
ture possibly undergoing a transition from monoclinic to 
triclinic symmetry PJJ . However, up to date the phys- 
ical mechanism underlying the rising of ferroelectricity 
is largely unknown. KhomskiiQ has suggested that 
ferroelectric (FE) polarization is caused by a combina- 
tion of site-centered and bond-centered charges between 
Fe 2+ /Fe 3+ ions in the CO B-site Fe4-tetrahedron pat- 
tern. His model, based on the structure of Ref.Q as- 
sumes each tetrahedron to show a "3:1" CO arrange- 
ment (three Fe 2+ and one Fe 3+ ions in a tetrahedron, 
or viceversa), at variance with Anderson's criterion, [l2| 
where each tetrahedron shows a "2:2" pattern (two Fe 2+ 
and two Fe 3+ ions). The 2:2 arrangement was found 
from density functional theory (DFT) [13] on the cubic 
Fd3m structure, distorted by X3 phonon mode. More- 
over, the 3:1 arrangement was recently obtained within 
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DFTpjJ Qj, LL5[ on a monoclinic P2/c structure. In ad- 
dition, a mixed CO pattern (25% of 2:2 and 75% of 3:1 
tetrahedra) was suggested to occur in a base-centered 
monoclinic Cc structure [16[. Our DFT results, in agree- 
ment with Ref.[16[, predict the Cc structure to be the 
ground state, consistently with recent results from reso- 
nant X-ray scattering [l7|. As discussed in Ref.[H, [~L9^ 
the Cc symmetry is stabilized by a delicate balance of 
different effects: Coulomb repulsion, entropy, [l2j , Fermi- 
surface nesting; [18] leading to a [001] charge-density wave 
(CDW),[l9[ oxygen breathing modes, etc. 

Recent experimental data showing real time FE switch- 
ing in magnetite epitaxial thin films were found to be in 
good agreement with our DFT calculations for the po- 
larization P in the Cc FE structure. [2l| In this letter, 
we focus on the proof of charge ordering as the micro- 
scopic origin of ferroelectricity in Fes04, by comparing 
two monoclinic P2/c and Cc lattice structures. Since 
the P2/c structure shows inversion symmetry, the total 
P must cancel out, so that at most antiferroelectricity 
may occur; on the other hand, the Cc structure is ferro- 
electrically active due to the lack of centrosymmetry. 

Electronic structure calculations were performed 
using the "Vienna Ab initio Simulation Package 
(VASP)",[23] where the projector-augmented-wave po- 
tentials, the generalized gradient approximation to the 
exchange-correlation potential [24| plus an effective on- 
site Coulomb interaction U were used.[25| In order to 
compare the total energies, the FE polarization and other 
relevant properties, we used - for both paraelectric (PE) 
P2/c and FE Cc states - the primitive cell of the base- 
centered monoclinic Cc lattice (with 112 atoms/cell) 
and experimental lattice parameters [20j. Since the dif- 
ference in experimental lattice constants between the 
PE and FE phases is ~0.1%, the same lattice param- 
eters in both unit-cells were used. The polarization vec- 
tor P is described with the conventional lattice vectors 
(a, 6, c). Internal atomic coordinates were optimized in 
both P2/c and Cc, starting from experimental Wyck- 
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FIG. 1: Ionic structure of Fe octahedral sites in P2/c (left) 
and Cc (right) cells. Orange and blue balls show Fe 2+ and 
Fe 3+ ions, respectively. Fe4 tetrahedra of 2:2 and 3:1 CO 
patterns are highlighted by yellow and green color planes, re- 
spectively. Electric dipole moments caused by charge shifts 
are indicated by red arrows. The lattice vectors in the prim- 
itive unit cell (x,y,z) and the conventional cell (a,&,c) are in- 
dicated. 



off parameters [3, [20j. The conventional atomic posi- 
tions of the Cc lattice were displaced by (-1/8, 0, 0) to 
fit into the P2/c structure, so as to have the pseudo- 
inversion center at (0, 0, 0). We focus on the P2/c and 
Cc structures since they are both proposed in diffraction 
experiments d, [20[ and show the lowest DFT total ener- 
gies compared to other proposed symmetries, i.e. Pmca 
and Pmc2i . jlq Note that the base-centered monoclinic 
lattice of Cc (a = b ^ c, a ^ (3 ^ 7 7^ 90°) is almost 
identical to a triclinic lattice (experimentally suggested 
as the FE state), the only difference being the a = b 
condition. A cutoff energy of 400 eV for plane waves, 
4x4x2 Monkhorst-Pack fc-point grid in the Brillouin 
zone, a threshold of atomic forces of 0.03 eV/A were used. 
An effective Coulomb energy U = 4.5 eV and an ex- 
change parameter J = 0.89 eV (as in Ref . [TBI] ) were used 
for Fe-d states, although some calculations with different 
values of U were performed (see below). Fe-3p 6 3d 6 4s 2 
and 0-2s 2 2p 4 electrons were treated as valence states. 
For Fe ions, the ferrimagnetic configuration was consid- 
ered, with all octahedral Fe sites as up-spin sites and all 
tetrahedral Fe sites as down-spin sites. Spin-orbit cou- 
pling was neglected. 

In terms of relevant structural and electronic proper- 
ties, our results are similar to the previous study by Jeng 
et. al. [Hj], so in this letter we will focus only on FE 
properties. As shown in FigJTJ octahedral Fe sites are 
located in xy planes with z = i/8 (i=0...7). The PE state 
has {E, C 2b + (0, 0, 1/2), J, a 2b + (0, 0, 1/2)} symmetries 
and the FE state has {E, a 2 b + (0,0,1/2)} symmetries 
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FIG. 2: Charge/orbital ordering of Fe minority tig states in 
PE configuration. Blue (orange) large balls show Fe 3+ (Fe 2+ ) 
ions. Red small balls show O ions. 



in a conventional base centered monoclinic cell so that 
there are two equivalent atoms (cfr. B12 and B12' sites 
in FigH]). Note that i) the mirror symmetry along with 
the translation vector forbids any net polarization along 
b and finite P is allowed only along the a and c directions; 
ii) the translation vector is relevant for the CDW stabi- 
lization and for opening the energy gap. [19] As discussed 
m 

Ref. mint 

the PE state shows entirely a 3:1 tetra- 
hedron CO arrangement whereas the FE state shows a 
mixed pattern. The difference between the two CO dis- 
tributions (see Fig. [2] for the P2/c and Fig. [3] for the 
Cc) can be understood when assuming a charge "shift" 
from B12 to B14 site and, in the upper part of the cell, 
from B12' to B14', all the other sites keeping their va- 
lence state. Each charge shift creates two 2:2 CO tetra- 
hedra, so as to form four 2:2 tetrahedra (cfr. FigH}. The 
resulting CO pattern lacks inversion symmetry, therefore 
allowing FE polarization. The CO rearrangement implies 
a change in the local breathing mode of O ions, which are 
driven away (attracted) by the substituted Fe 2+ ( Fe 3+ ) 
ion at B12 (B14) site, with a displacement of ~ 0.1A(cfr 
small black arrows in Fig. [3]). Similarly, upon charge 
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TABLE I: P along a, b and c axes (in /iC/cm 2 ), calculated 
with different approaches (see text). 

Pfierry P?CM Prfip 

(-4.41, 0, 4.12) ( -4.20, 0, 5.27) (-4.05, 0, 5.73) 



model (Pdip) where we summed up the two dipole mo- 
ments located at sites where charge-shifting occurs (i.e. 
from B12 to B14 site) with nominally one electron (see 
red arrows in Fig. [I]). The consistency of these values, 
shown in Table H implies two important facts: (i) P is 
induced largely by the CO rearrangement (i.e. charge 
shifts at few sites) but not much by the ionic displace- 
ments, as evident by the similarity between Ppcm and 
Pdipi (ii) since the Berry phase value is equivalent to 
the sum of Wannier- function (WF) centers [26], the sim- 
ilarity between P Berry and ~Ppcm suggests the WF of 
Fe-d electron to be mostly centered at the ionic sites, no 
matter how the WF is delocalized and hybridized with 
surrounding O-p orbit als. 



TABLE II: Charge separation (cs, i.e. difference of <i-charges 
between Fe 2+ and Fe 3+ ions in the atomic sphere with lA 
radius) and the corresponding FE P Berry (in /iC/cm 2 ) vs 
Coulomb repulsion U (J is fixed to 0.89 eV). 
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(-4.41, 0, 4.12) 


(-4.42, 0, 4.81) 


(-4.33, 0, 5.07) 



FIG. 3: Charge/orbital ordering of Fe minority tig states in 
FE configuration (notation as in Fig[2j . Black arrows show 
atomic displacements with respect to the PE reference state. 



shifts of the Fe ti g electron between sites, some of the 
Fe 3+ ions (at B9 and Bll sites) move towards the Fe 2+ 
ions, maybe reminiscent of Khomskii's mixed bond-/site- 
centered charge mechanism. [?} However, as far as fer- 
roelectricity is concerned, when focusing on a Fe chain 
along [100], these movements average out so as not to 
give a net contribution to P. 

As for calculated P (cfr Table [Tj), the Berry phase 
approach [2^ predicts quite a large polarization (at least 
compared to other improper multiferroicsp],[2[), its direc- 
tion lying in the ac mirror plane. As noted previously [2lj. 
the DFT results are in excellent agreement with recently 
reported experimental values for magnetite thin films 
(P ~ 5.5 fiC/cm 2 in the ab plane with the c compo- 
nent not measured) as well as with earlier experiments 
on single crystals [lOj: P a = 4.8 11C / cm 2 and P c = 1.5 
liC/cm 2 . 

To deepen our analysis, we compare the value of 
P Berry with the value estimated from a point charge 
model, (Ppcm), assuming the nominal valence for every 
ion, i.e. "full" charge disproportionation, i.e. 2+ and 3+ 
for charge-ordered Fe, 2- for O) as well as with a simple 



Our picture based on the WF centers is also con- 
firmed by the results of P Berry values upon varying the 
U Coulomb parameter. As shown in table [III upon in- 
creasing U and keeping the atomic configuration fixed to 
that obtained for U = 4.5 eV, the charge separation be- 
tween Fe 2+ and Fe 3+ is enhanced. Note that the value of 
P does not change rapidly with U, suggesting the centers 
of the WF not to move significantly far from the Fe 2+ for 
all U values. Also, in the limit of extemely large J7, we 
expect a "full" charge disproportionation to occur, caus- 
ing P Berry to become progressively closer to Ppcm (as 
confirmed by Table HI|) . 

For further insights, an "adiabatic" path is set up to 
connect the PE and FE state by displacing all ions lin- 
early with a scaling parameter A (i.e. A=l for full FE dis- 
placements and A=0 for the initial PE structure). The 
FE structure with opposite P was built starting from 
the PE structure and with A=-l, i.e. with displacements 
opposite to the structure considered so far. By anal- 
ogy with the above discussion for positive A values, two 
charge shifts between B10 and B6 (B10' and B6') occur 
as a transition from the PE to the "negative" FE phase. 
The charge separation (cs) is calculated by comparing the 
"muffin-tin" charge on the sites where the shift occurs in 
going from PE to FE (i.e. B12 - B14 sites when A > 
0). The results are reported in Fig HI The total energy 
shows a double valley structure, typical for ferroelectric- 
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ity, with a deep global minimum at the FE state, along 
with a very shallow local minimum at the PE state. In 
going from the PE to the FE state, we show that the 
charge shift rather suddenly occurs around A=0.5, when 
P is strongly enhanced and the total energy reduced (cfr. 
FigHla) and b)). Moreover, the origin of ferroelectricity 
as driven by charge rearrangements on different sites is 
clearly confirmed by the trend of P x closely following that 
of cs along the path ((cfr. FigJUc)). 

In summary, we have shed light on the microscopic ori- 
gin of ferroelectric polarization in insulating magnetite 
by analyzing the differences between two charge-ordered 
states: the P2/c PE state and the Cc FE state. Ferro- 
electricity is induced by a non-centrosymmetric charge- 
ordering and the polarization is primarily caused by local 
dipoles at selected octahedral sites, pointing to a pic- 
ture of ferroelectricity mostly based on localized "charge 
shifts". Our calculations show that FesC^ can be consid- 
ered as a prototypical case of charge-order driving multi- 
ferroicity with relatively large values for the electric po- 
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